Superfluid Interfaces in Quantum Solids 
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One scenario for the non-classical moment of inertia of solid *He discovered by Kim and Chan 
[Nature, 427, 225 (2004)] is the superfluidity of micro-crystallite interfaces. On the basis of the most 
simple model of a quantum crystal — the checkerboard lattice solid — we show that the superfluidity 
of interfaces between solid domains can exist in a wide range of parameters. At strong enough inter- 
particle interaction, a superfluid interface becomes an insulator via a quantum phase transition. 
Under the conditions of particle-hole symmetry, the transition is of the standard U{1) universality 
class in 3D, while in 2D the onset of superfluidity is accompanied by the interface roughening, driven 
by fractionally charged topological excitations. 
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Recent observation by Kim and Chan of non- 
classical moment of inertia (NCMI) of "^He at pres- 
sures significantly higher than the sohdification point 
is a breathtaking result, especiaUy striking in view 
of the theorem-like theoretical arguments against ex- 
istence of commensurate supersolids :2i. The fact of 
commensurability — equivalently, one may put it as the 
fact of absence of vacancies, or interstitials, or both — of 
the equihbrium sohd "^He at T = is supported by an ex- 
tensive experimental work over the past several decades 
(for review, see, e.g., 0); as well as by the most recent ex- 
perimental and numeric studies [J, |a| ■ The commensura- 
bility of solid ''He rules out NCMI based on Bose-Einstein 
condensation of vacancies 6] . Two of us have proposed 
recently that NCMI might be due to the superfluidity 
of interfaces between ''He crystallites. At present, the 
weak point of this hypothesis is the absence of a theoret- 
ical analysis and/or direct experimental evidence of the 
superfluidity in the walls separating insulating domains. 

The problem of interface superfluidity in a quantum 
solid is of signiflcant general interest on its own, being 
potentially relevant not only to the solid '''He polycrystal, 
but also to the properties of domain walls in spin arrays 
and ultracold atoms in optical lattices. 

In this Letter, we present a proof-of-principle study 
of superfluidity in interfaces between insulating domains 
with broken translation symmetry (solids). We address 
the problem by studying the checkerboard lattice solid 
(CB). We start with giving a simple illustrative theoret- 
ical argument that at least under certain limiting condi- 
tions the domain wall in our system has to be superfluid. 
Our numeric simulations of 2D and 3D models reveal su- 
perfluidity of the CB domain walls in a large range of 
parameters. We pay special attention to the study of the 
superfluid (SF)-insulator (I) quantum phase transition 
in the interface. In 3D, the transition turns out to be in 
the U{1) universality class. In 2D, we conclude that the 
I-to-SF transition in the wall is driven by proliferation of 
topological excitations that carry the fractional particle 
charge 1/2, as well as the quantum of the interface shift in 



the transversal direction. The latter circumstance results 
in an interesting effect: The transition is accompanied by 
the wall roughening. 

The simplest system featuring both superfluid and 
CB phases is that of the hard-core lattice bosons with 
nearest-neighbor repulsion, at half-integer filling factor 
(see, e.g., \f] and references therein). The model can 
be exactly mapped onto spin-1/2 XXZ antiferromag- 
net, which leads to the following correspondence. The 
CB phase is equivalent to the easy-axis antiferromagnet 
(characterized by the broken Z2 symmetry), while the 
SF phase is identified with the easy-plane antiferromag- 
net (characterized by broken U{1) symmetry). Corre- 
spondingly, in terms of the Neel vector, S, the CB order 
parameter is M = Sz, while the SF order parameter is 

= Sx + iSy. Generically, the groundstate of the model 
is either SF or CB, depending on the Hamiltonian param- 
eters. There is also a special S'C/(2)-symmetric Heisen- 
berg point. In 2D and 3D groundstates of the Heisenberg 
Hamiltonian the SU (2) symmetry is broken, so that the 
vector S is non-zero and can point at any direction. 

Let us take now an easy-axis Hamiltonian that is very 
close to the Heisenberg point (T = 0), and create two 
large domains, Sz = M and Sz — —M. What is the 
structure of the domain wall? Being close to the Heisen- 
berg point, we are forced to conclude that the wall is 
very thick in the transverse direction, with the vector S 
well defied locally inside the wall and evolving smoothly 
from (0, 0, M) to (0, 0, ~M) across the wall. The ener- 
getic cost is controlled by the closeness to the Heisenberg 
point and can be rendered arbitrarily small. In the mid- 
dle of the wall, S = {Sx, Sy,0), which means that the 
wall is characterized by broken U{1) symmetry implying 
superfiuidity in the bosonic case. As the system is taken 
deeper into the solid state in the bulk, it becomes energet- 
ically favorable to suppress the module of S to zero in the 
wall, which means an insulating state of the wall. Similar 
transformation of the domain wall — from the Bloch-type, 
where \S\ is finite everywhere, to the Ising-type, where 
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\S\ = in the middle of the wall — takesplace in classical 
ferromagnets close to the Curie point j^]- 

Once the superfluidity of a domain wall is established 
in the limiting case, one may expect that it can take place 
under more general conditions, especially in view of the 
following, almost obvious, energetic argument. Due to 
geometrical frustration, the energy cost to translate a 
particle along the wall is less than in the bulk. For ex- 
ample, in the case of the 2D CB solid, the particle jump 
increases the energy by ~ 3^, where V > is the in- 
teraction energy between two close neighbors. A particle 
at the wall [see, e.g., Fig.0Ja)] already has another par- 
ticle as its neighbor. After jumping one step along the 
wall, it acquires only two close neighbors, and the energy 
increase due to such a jump is only V. 

A point of concern, however, is that generically the 
transition from SF to CB is of the first order so that 
for a given system or range of parameters it may turn 
out that the interface between two insulating domains is 
always in the insulating state. 

To get an idea of how likely it is to get the interface 
between the two CB domains superfluid, we simulate 
a domain wall in the bond-current model — a discrete- 
imaginary-time analog of a {d+ l)-dimensional worldline 
representation of a quantum bosonic or spin system in 
d spatial dimensions 9]. The Hamiltonian of the model 
reads 

^ = ^EE'^n." + f E J^,rJm,r. (1) 
n a=l (n,m) 

Here the integer vector n — (rii, . . . , Ud, rir) labels sites 
of the {d+ l)-dimensional cubic lattice, a = 1, . . . , d enu- 
merates the spatial directions, and r denotes the tempo- 
ral direction. Currents Jn,a and J^t arc integers asso- 
ciated with bonds adjacent to the site n in directions a 
and T, respectively. The summation in the second term 
runs over all pairs of temporal bonds having a common 
plaquette. The configurations of bond currents are sub- 
ject to the zero-divergence constraint. Without loss of 
universality, we restrict the values of bond currents to 
Jn,a = 0,±1 and Jn,T = ±1- The zero-divergence con- 
straint then reads 

( J„,a + J„,-a) + ^ ( J„,r + Jn,-r) = , (2) 

a 

where the negative sign means the opposite direction, so 
that Jn,-a = — t/(n-a),Q ^nd Jn,-T = —J(ii-f)y, the hats 
stand for unit vectors in the corresponding directions. 

We simulate the model using Worm algorithm 
|lfl| |. in the range of parameters where the bulk is deep 
into the solid regime. To automatically create a do- 
main wall, we take a lattice with periodic boundary con- 
ditions and an odd number of sites in the spatial di- 
rection X (system sizes in all the other directions are 
even). The domain wall superfluidity manifests itself 
as a non-zero mean square of winding numbers ,11.] in 



the direction(s) parallel to the wall: {W^) 7^ 0, while 
(W^) = 0. More specifically, for a 3D system with the 
domain wall in the yz plane, the superfluid stiffness is 
given by ps = {Wy + W^)/2Lr, and likewise for the com- 
pressibility. In a 2D system with the domain wall in 
the y direction the statistics of winding numbers is es- 
sentially discrete, P(W^y) oc exp [—{Ly/Lr)Wy/2ps], and 
an appropriate estimator for the superfluid stiffness is 
p-^ = -2{Lr/Ly)\n[P{Wy = l)/P{Wy = 0)], and like- 
wise for the compressibility. Here and Ly are the lin- 
ear system sizes in the corresponding directions. In our 
simulations, we set Ly = = = L and = L + 1. 

In 3D, we found that for t < 1 the interfaces are 
never superfluid. Simulations at t — 1.3 revealed a first- 
order SF-CB transition in the bulk at p w 0.2 with the 
interface remaining superfluid till p becomes equal to 
Pc = 0.27115(5). In 2D, the simulations were performed 
at t = 0.9. The flrst-order bulk SF-CB transition was 
found at p « 0.5, while the interface becomes insulating 
only at Pc = 0.7633(5). 

A superfluid interface embedded into a d-dimensional 
solid is a (d — l)-dimensional superfluid, and the sce- 
nario for the quantum phase transition from SF to I 
state in such a system is interesting on its own. Our 
model, Eq. has a particle-hole symmetry. Thus, com- 
mensurability should play a key part in the criticality 
|12|. Fundamentally, there are two qualitatively differ- 
ent cases, depending on whether the interface is smooth 
or rough. If the interface is centered at a; = and x{p) 
is its instantaneous shape (p is the vector in the hyper- 
plane perpendicular to the x axis), then, by deflnition, 
{x{pf' ) ^ 1 for a smooth interface, while for a rough in- 
terface {x{pf') is macroscopically large (scales as some 
power of the system size). For a smooth interface, the 
CB environment plays a role of a periodic external po- 
tential that doubles the interface unit cell. This means 
that a smooth interface can be treated as a commensu- 
rate system with an integer (unity) filling factor. Its SF-I 
transition then corresponds to the superfluid-Mott insu- 
lator (MI) transition at integer filling, known to be of 
the U{1) universality class [Ijl- If the interface is rough, 
then the effect of the solid environment is averaged out 
by the zero-point fiuctuations of x[p) and the effective 
filling factor for the interface remains half-integer, with 
corresponding implications for the universality class of 
the SF-I quantum phase transition. 

In our simulations, we observe the smooth-interface 
scenario in 3D and the rough-interface scenario in 2D. 

In Fig. n we present the 3D data in the vicinity of the 
critical point, Pc — 0.2711(5). A very good data collapse 
with the critical exponent v k, 0.671 of the U{1) univer- 
sality class is indicative of the standard SF-MI scenario. 

In 2D, the interface forms a ID Luttinger liquid, the 
main characteristic of which is the dimensionless param- 
eter g = tt^/pJk, where pg and k are the ID superfluid 
stiffness and compressibility, respectively. In a Luttinger 
liquid with a flUing factor \/m (m is an integer) the 
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FIG. 1: Criticality for the interface in a 3D system. The data FIG. 2: Luttinger hquid parameter y as a function of p. Dot- 
for ps as a function of the parameter p is consistent with the ted hues are to guide the eye. The solid line is the KT ex- 
exponent u ^ 0.671 of the 17(1) universality class. trapolation to the infinite system size. 



SF-I (Kosterlitz-Thouless type) transition takes place at 
g = gc — 2/m~ 14]. As we argued above, a smooth in- 
terface implies m = 1 and, correspondingly, criticality at 
gc = 2. The results of our simulation show that the SF-I 
transition actually takes place at gc = 1/2, correspond- 
ing to the half-filling case, implying the rough interface 
scenario. 

As is always the case with SF-I transitions in ID 
systems, a brute-force numeric observation of the crit- 
ical gc is problematic in view of the exponentially diver- 
gent correlation length. We thus need to perform the 
finite-size analysis of the data using Kosterlitz-Thouless 
renormalization-group flow: 

r9{L2)/ga 

-5-— — — =41n(Li/L2). (3) 

Here g{L) is the Luttinger parameter g as a function 
of the system size, ^ is an i-independent microscopic 
parameter (which is an analytic function of p). At a 
given p, the value of ^ is obtained with Eq. (j^J from 
numeric values of g{Li) and g{L2). The consistency 
with the Kosterlitz-Thouless renormalization-group flow 
is checked by the data collapse for different pairs of sys- 
tem sizes and also by the shape of the curve ^{p) which 
should look as a straight line in the vicinity of the criti- 
cal point, in contrast to the g{L,p) curves — as functions 
of p, at large enough L. These curves should demon- 
strate a considerable curvature consistent with a (slow) 
evolution, as i — > oo, towards the jump at pc from 
g = 1/2 to g — Q. These features are seen in Figs. 
|3| After extracting the function £,{p), the macroscopic 
g{p) = g{L — > oo,p) limit is obtained — in accordance 
with Eq. ©—from (g/gc) (Ing/gc - = -1- 

The critical value g — 1/2 implies a rough interface, 
since roughening is apparently the only mechanism of 
eliminating the effect of broken translation symmetry in 
the bulk. There is also a strong argument in favor of si- 
multaneous appearance of superfluidity and roughening 
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FIG. 3: Parameter ^ for different data sets as a function of 
p. The solid line is a linear fit. The errors are of the order of 
the symbol size. 



in the ID interface. There is little doubt that deep in the 
insulating phase the interface becomes smooth (T = 0), 
since roughening costs finite potential energy that domi- 
nates over the kinetic energy in this limit. In our simula- 
tions, we see this as the effect of "freezing" of the interface 
position at large enough p. [By its nature, a rough inter- 
face experiences local fluctuations that gradually lead to 
its global drift.] The zero-point roughening fluctuations 
in the smooth phase are due to the specific solitons il- 
lustrated in Fig. 21 These solitons shift the position of 
the interface in the x direction by one step. It is also 
seen that they carry a topological charge associated with 
shifting by one lattice period the checkerboard density 
wave along the interface (y direction). These quasiparti- 
cles also carry a particle charge ±1/2, as it follows, e.g., 
from the fact that a single-particle hopping event trans- 
lates the soliton by two lattice spacings, see Fig. 0] In 
view of their fractional particle charge it is conventional 
to call these solitons spinous 15]. It turns out that in the 
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FIG. 4: A sketch of the wall roughening mechanism due to spinous, (a): a "frozen" wall (the arrow shows a single-particle 
hopping event, which generates the configuration (b), and the dashed lines are to facilitate the wall gazing); (b): a configuration, 
featuring a pair of spinous; (c): with a single-particle hopping over one lattice site the spinon is shifted by two lattice periods, 
which means that the spinon particle charge is 1/2. 



insulating phase spinons are the lowest particle-charge 
carrying elementary excitations. — The snapshots of the 
wordline configuration cross-sections in planes perpen- 
dicular to the T direction reveal structures identical to 
those of Fig. 0| It is reasonable to assume then that 
the transition from insulating to supcrfluid phase is due 
to the proliferation of spinons, in analogy to the super- 
fluid transition in a standard ID checkerboard solid, also 
driven by spinons. Hence, we arrive at a picture where 
both superfluid and roughening transitions occur simul- 
taneously being driven by proliferation (condensation) of 
one and the same quasiparticle mode. 

In conclusion, we have demonstrated that an interface 
layer in a normal solid may exhibit superfluidity in a wide 
range of parameters. This result may be of direct rele- 
vance to NCMI of solid ''He discovered by Kim and Chan 
0, supporting the interpretation in terms of the super- 
fluidity of micro-crystallite interfaces. We have studied 
numerically superfluid-insulator quantum phase transi- 



tions in particle-hole symmetric interfaces in 2D and 3D 
models of the lattice checkerboard solid. In 3D, the tran- 
sition is in the U{1) universality class implying that the 
interface is smooth. In 2D, where the interface is a ID 
Luttinger liquid, we observe a Kosterlitz-Thouless type 
transition at Luttinger-liquid parameter g — 1/2, which 
implies that the interface is rough on the superfluid side. 
We argue that the ID interface becomes smooth simulta- 
neously with becoming insulating, since the onset of su- 
perfluidity and roughening are due to proliferation of the 
same quasiparticles which (i) have particle charge ±1/2, 
(ii) represent defects in the checkerboard order, and (iii) 
are kinks shifting the interface in the perpendicular di- 
rection by one lattice spacing. 
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